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This is a pedagogical review of the subject of linear polymers on deterministic finitely 
ramified fractals. For these, one can determine the critical properties exactly by real-space 
renormalization group technique. We show how this is used to determine the critical 
exponents of self-avoiding walks on different fractals. The behavior of critical exponents 
for the n-simplex lattice in the limit of large n is determined. We study self-avoiding 
walks when the fractal dimension of the underlying lattice is just below 2. We then 
consider the case of linear polymers with attractive interactions, which on some fractals 
leads to a collapse transition. The fractals also provide a setting where the adsorption of 
a linear chain near on attractive substrate surface and the zipping-unzipping transition of 
two mutually interacting chains can be studied analytically. We also discuss briefly the 
critical properties of branched polymers on fractals. 

1. INTRODUCTION 

The problem of self-avoiding walks ( SAWs) on lattices provides perhaps the simplest 
geometrical model of equilibrium critical phenomena ( i.e., non-trivial power-laws in the 
behavior of different quantities in a system in thermal equilibrium). The two other familiar 
examples of geometrical models showing phase transitions, the percolation problem and 
a system of hard particles (spheres or rods), both involve more complex geometrical 
structures. The model of SAWs captures the important macroscopic features of polymers 
in solution, and is closely related to other models of phase-transitions in statistical physics 
like the Ising model [|T], and can also be seen as the n —>■ limit of the n- vector model [ 
121. Given the many technological applications of polymers, and the importance of SAWs 
as a model of critical phenomena, it is not surprizing that the model has attracted a lot 
of attention in the last sixty years. Several good reviews are available summarizing our 
current understanding of this problem [E]. 

The SAW problem is clearly trivial in one dimension. In spite of the large number of 
papers related to this problem, an exact solution of the problem has not been possible so 
far, for any non-trivial case. In two dimensions, the exact value of the growth constant 
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is known for the hexagonal lattice [E], and has been conjectured for the square lattice 
from large-order exact series expansions [IH]. The exact value of the critical exponents 
in two dimensions is known from conformal field theory [IB]. Even for dimensions d > 4, 
where the critical exponents are known to take mean-field values [Ej, a full solution has 
not been possible. 

Given the paucity of exact solutions in this area, it seems reasonable to look for some 
artificially constructed graphs, e.g. fractals, for which an exact solution can be found. 
This solution then, can be considered as an approximation to the original problem. The 
advantage of this approach, over other ad-hoc approximations like the Flory approxima- 
tion, is that one is assured of well-behavior requirements like the convexity of the free 
energy, and avoids problems like getting two different values for a quantity ( e.g. the 
pressure for hard-sphere systems), if one calculates it in two different ways within the 
same approximation. 

Another motivation for the study of SAWs on fractals comes from the need to under- 
stand how the critical exponents depend on the dimension and the topological structure 
of underlying space. In the formal techniques like e-expansions, one treats the dimension 
of space (i as a parameter that can be varied continuously, but such techniques involve 
a formal analytical continuation of various perturbative field theory integrals of the type 
/ d'^k.. to all real values of c? [|H]. However, this formal approach does not give any simple 
answer to the question, "What is the Ising model in 3.99 dimensions ?". On the other 
hand, fractals are explicitly constructed spaces of non-integer dimension, and one can 
construct fractals with dimension close to any given real number d. One can thus study 
how critical exponents change by changing the geometry of the underlying space. Inter- 
estingly, the results from the formal e-expansion techniques do not match with those from 
explicitly constructed fractals. The main reason seems to be the assumption of transla- 
tional invariance in the former, while the explicitly constructed fractals do not have this 
property. In fact, this assumption seems to be problematic for spaces with non- integer 
dimensions, and leads to pathologies such as non-positive measure For a longer 

discussion of these issues, see [ ITU] . 

A third reason for the study of SAWs on fractals is that these provide excellent peda- 
gogical examples of application of renormalization group techniques to the determination 
of critical exponents. If one tries to implement the renormalization transformation to 
some system like the Ising model in two dimensions, one eventually generates an infinity 
of additional multi-spin or longer range couplings, which are presumably irrelevant, and 
are neglected to get renormalization equations in terms of a finite number of variables. 
The justification given for doing this usually involves too much 'handwaving' for a dis- 
cerning student. In case of fractals, one gets exact renormalization equations in terms of 
a finite number of coupling constants. One can then study these in detail, and explicitly 
work out their trivial and non-trivial fixed points, stable and unstable directions, basins 
of attractions of fixed points, critical exponents in terms of eigenvalues of the linearized 
transformation etc., and learn the use of renormalization group techniques for determining 
critical properties. 

Finally, one can treat the fractal lattice as a simple model for a disordered substrate on 
which the polymer is adsorbed, and use the exact results found for polymers on fractals 
to develop some understanding about real experimental systems. But for this, this article 
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would not have been included in this volume. 

The simplest of graphs for which the SAW problem is analytically tractable is the Bethe 
lattice, for which the exact solution is trivial, as the graph has no loops. Only slightly 
more complicated is the case of Husimi cactus graphs, which are like the Bethe lattice on 
the large scale, but have small loops [Fig. These do change some properties like the 
connectivity constant for the walks, but the large-scale properties of walks on such lattices 
are the same as on the Bethe lattice. In this chapter, we will discuss fractal graphs which 
have loops of arbitrarily large size, and for which the critical exponents are different from 
the mean-field values. 

This chapter is organized as follows: In section 2, we start with the precise definitions 
of the fractals we study. We also define the SAW generating functions, their annealed 
and quenched averages and the critical exponents of self-avoiding walks. In section 3, we 
discuss the renormalization group treatment of the SAWs on the 3-simplex lattice in detail. 
We explicitly construct the recursion equations, and show how the critical exponents 
of SAWs can be determined from an analysis of the linearized renormalization group 
transformation near the fixed point. The dense phase of the polymer, and the fluctuations 
of the number of rooted walks are also discussed briefly. In section 4, we show how 
this treatment is generalized to other fractals where one has more complicated recursion 
equations than the simplest case of 3-simplex lattice.. We also discuss the behavior of 
critical exponents on the ra-simplex in the limit of large n. In section 5, we discuss 
the SAWs on the Given-Mandelbrot family of fractals, and study the behavior of critical 
exponents of walks when the fractal dimension of the lattice tends to 2 from below. We use 
finite-size scaling theory to determine how the structure of the renormalization equations 
depends on the parameter b defining the fractals. We develop a perturbative expansion 
for the critical exponents valid for large b when the fractal dimension of the lattice is just 
below 2. In section 6, we discuss the collapse transition of linear polymers with attractive 
self-interaction, and the tricritical ^-point. We also discuss other intermediate 'semi- 
compact' phases that is seen on some fractals. In section 7, we show that much of our 
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Figure 2. The recursive construction of the n-simplex fractal for n = 5. (a) The first 
order graph (b) the graph of a (r + 1) order triangle, formed by joining n r-th order 
graphs, shown as shaded blobs here (c) the second order graph. 



treatment of linear polymers can be extended to branched polymers. One can determine 
the critical exponents 6 for average number of branched polymer configurations per site, 
and u for the mean polymer size using the RSRG techniques. Interestingly, one finds 
that for GMh fractals with b > 2, the average number of branched polymers per site 
increases as exp(a?T, + bn^), with < < 1, and the leading form of correction to the 
exponential growth is not a power-law correction. This stretched exponential form arises 
from the presence of favorable and unfavorable regions on the fractal lattice. In section 
8, we discuss a polymer with self-attraction near an attractive surface. In this case, 
there is a competition between the collapse transition in the bulk of the fractal, and the 
tendency of the polymer to stick to the surface. We discuss the qualitative features of 
the phase-diagram, and critical exponents. The phase-diagram is somewhat complicated 
when the next layer interaction is included, and provides a good pedagogical example for 
the study of higher order multi-critical points. One can study even more complicated 
systems. Section 9 contains a discussion of two mutually interacting linear chains, as a 
model of the zipping -unzipping transition in the double-stranded DNA molecules. Section 
10 contains some concluding remarks, and we mention some open problems which deserve 
further study. 

2. DEFINITIONS 

We start by giving definitions of the family of fractals that we study in this article. All 
these fractals have a finite ramification number, and are defined recursively. 

The first family we shall discuss is the n-simplex family, defined for all positive integers 
n > 2. The n = 3 case was first defined by Nelson and Fisher [ TT], who called it 
truncated tetrahedron lattice. The construction was generalized for arbitrary n in [I12j. 
The recursive construction of the graph of the fractal is shown in Fig. |^ The first order 
graph is a single vertex with n bonds. In general, the r-th order graph will have rf~^ 
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Figure 3. The recursive construction of the Given-Mandelbrot fractal for 6 = 3. (a) The 
graph of first order triangle, (b) the graph of a (r + 1) order triangle, formed by joining 
b{b + l)/2 r-th order triangles shown as shaded triangles here (c) graph of the 2nd order 
triangle. 



vertices, and {rf + n)/2 bonds. Of these {rf — n)/2 are internal bonds, and there are 
n external bonds, that are used to connect this graph to other vertices to form bigger 
graphs. To form a graph of (r + l)-th order, we take n graphs of r-th order, and join 
them to each other by connecting a dangling bond of each to a dangling bond of the other 
r-order subgraphs. There is one dangling bond left in each of the subgraphs, and n bonds 
altogether. These form the dangling bonds of the (r -|- l)-th order graph. We let r tend 
to infinity to get an infinite graph. The fractal dimension D„ of this graph is easily seen 
to be log n/ log 2. It was shown in [ 12j that the spectral dimension (i„ of this graph is 
21og?T,/ log(?T, -I- 2). 

The second family we shall consider is the Given-Mandelbrot family of fractals [ I13j. 
defined in Fig. El Members of the family are characterized by an integer b, with 2 < b < 
oo. We start with a graph with three vertices and three edges forming a triangle. This is 
called the first order triangle. To construct the graph of the (r + l)-th order triangle, we 
join together graphs of b{b + l)/2 triangles of r-th order, (i.e. identify corner vertices) as 
shown in the figure, to form an equilateral triangle with base which is b times longer. We 
shall call the resulting graph the GMb fractal. 

It is easy to see that the number of edges in the graph of the r-th order triangle is 
3[b{b+ l)/2Y-\ The fractal dimension is A = log[6(6 + l)/2]/log6. For b = 2,3,4.., 
these values are 1.5849,1.6309,1.6609... respectively. The spectral dimension of the 
graph can also be calculated exactly for general b [[H|- The values of (ib for 6 = 2 to 10 
are listed in [CHI. For large b, db tends to 2, and the leading correction to its limiting 
value is given hj db ^ 2 — log log 6/ log 6 [ITB]. 

In fact, following Hilfer and Blumen [|^, one can define a general fractal family of 
Sierpinski -like fractals, to be called HB{b,d) family here, such that the the n-simplex 
corresponds to HB{2,n — 1), and the GMb fractal corresponds to HB{b,2). Here the 
basic unit is a (d-fl) -simplex graph, and one makes the {r + l)-th order graph by making 
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a simplex of b layers of the r-th order graphs. The construction is illustrated in Fig. I^for 
the fractal HB{3,3). 




Figure 4. First two stages in the iterative construction of the Hilfer-Blumen fractal 
HB{3,3). This has a fractal dimension d = Iogl0/log3. 



The third family of fractals we shall use is a generalization of the modified rectangular 
lattice [1^. Here the first order graph is a set of 2'^ vertices forming a (i-dimensional 
hypercube. In general, the r-th order graph has 2^"'"'^"^ vertices in the shape of a cuboid, 
with each corner vertex of the cuboid having an extra dangling bond to connect to outside. 
The graph of an (r + l)-th order cuboid is formed by taking two r-th order cuboids, and 
joining the opposite faces by 2°'~^ bonds. The direction of the face that is selected for 
this is changed cyclically at the next order, so that the lengths of different sides of the 
cuboid are within a factor 2 of each other. It is easily seen that the fractal dimension 
of this graph is d, and it can be shown that the spectral dimension is 2(1 — 2"*^). The 
construction is illustrated in Fig. El for the case d = 2. 

The determination of the generating function for the linear polymers on these fractals 
follows the treatment of [El^]. Let Pn{N) be the number of distinct simple polygons of 
perimeter n on a graph with N total number of sites in the graph, different translations 
of the polymer being counted as distinct. [The use of same symbol n for the n-simplex 
and n-stepped walk should not cause any confusion, as it is clear from context which is 
intended.] For large A^, this number increases linearly with A^. We then define P„ as the 
average number of polygons of perimeter n per site by 

Pn = LimN-^oo^-^j^; (1) 
and 

oo 

P(X) = Y: ^nX", (2) 
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Figure 5. The recursive construction of the modified rectangular lattice for d = 2. (a) 
The graph of first order square. (b) the (r + l)-th order graph, formed by joining 2 r-th 
order graphs shown as shaded rectangles here (c) graph of the 5th order rectangle. 



For large n, P„ varies as fi"'n°'~^, where is a constant, and a is a critical exponent for 
the walks. This corresponds to P{x) ~ (1 — ^x)'^''^, for x tending to from below. 
Note that P{x) is a finite-degree polynomial in x for the graphs like the Husimi cactus 
(Fig.P). 

We similarly define C„, the average number of open walks of length n, as the average 
over all positions of the root, of the number of n-stepped self-avoiding walks. For large n, 
this varies as /i"'n'^~^, which defines the critical exponent 7. The corresponding generating 
function 

00 

C{x) = J2 CnX"" (3) 

n=l 

varies as (1 —xfi)~'^ as x tends to from below. We also define the critical exponent p 
for SAWs by the relation that the average end-to-end distance for walks of length n varies 
as for large n. 

Here in the context of disordered systems, we imagine that the polymer can freely 
move over all space, and the averages calculated correspond to annealed averages, as we 
are averaging the partition function of SAWs over different positions of the root. The 
logarithm of the number of polymer configurations is the entropy, and one can also define 
the equivalent of quenched averages, where one averages, not the numbers of walks, but 
logarithms of numbers of rooted walks, over different positions of the root. These are more 
difficult to determine. We will indicate how these can be handled in sec. 3. 

3. RENORMALIZATION EQUATIONS FOR THE 3-SIMPLEX FRACTAL 

The procedure to determine the critical behavior of SAWs on fractals is simplest to 
illustrate using the 3-simplex as an example. We will keep the presentation informal 
(hopefully not inaccurate). Readers who prefer a more formal approach, may consult [ 
ITU] . We discuss the calculation of annealed averages first. 

3.1. Calculation of annealed averages 

Consider one r-th order triangular subgraph of the infinite order graph. It is connected 
to the rest of the lattice by only three bonds. Our aim to sum over different configurations 
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of the SAW on this subgraph, with a weight x for each step of the walk. These can be 
divided into four classes, as shown in Fig. IHl Here A^^^ is the sum over all configurations 
of the walk within the r-th order triangle, that enters the triangle from a specified corner, 
and with one endpoint inside the triangle. B^"^^ is the sum over all configurations of walk 
within the triangle that enters and leaves the triangle from specified corners. In C^^\ 
the walk enters and leaves the subgraph from specified vertices, and reenters afterwards 
from the third corner, and has one endpoint inside the triangle. D^^'' consists of sum over 
configurations that have both endpoints of the walk inside the triangle, but part of the 
walk is outside the triangle. We can write down the values of these for r = 1 immediately. 

v4(^) = v^, 5(1) = X, C(i) = = 0. (4) 



We define the order of a closed polygon on the infinite 3-simplex graph to be r, if the 
polygon is completely contained inside an r-th order triangle, but not inside any (r — l)-th 
order triangle. It is easy to see that sum over all r-th order polygons within one such 
triangle is B^"^^^^ . The number of sites in the (r + l)-th order triangle is S*", hence we get 

P(x) = ^3-"fi(^) . (5) 

r=l 

The sum over open walks can be expressed similarly 
C(x) = S^lSA^") + 35("M(") + W^'^ (6) 

r=l 

It is straight forward to write down the recursion equations for these weig 
C*-'"'*"^'' and Z}(^+^) in terms of the values at order r. For example. Fig. [7| shows the only 
two possible ways one can construct a polymer configuration of type B. It is easy to verify 



that the resulting equations are [ fTTj 

A^'-+^^ = A{l + 2B + 2B^) + C{2B'^)] (7) 

Qir+l) ^^2^ ^3. (g) 

= A52 + C(352); (9) 

jj{r+i) ^ (^2 ^ 2A'^B + AABC + QBC"^) + D{2B + SS^); (10) 



where we have dropped the superscripts (r) in the right-hand side of the equations to 
simplify notation. 

It is straightforward to determine the analytical behavior of the generating functions 
P{x) and C{x) from the equations (4)-(10). 

We note that the equation for i?('"+^) depends only on B^'^\ and does not involve the 
other variables. The recursion equation (jHl) has two trivial fixed points B* = 0, and 
B* = oo, and a non-trivial fixed point B* given by B* = (a/5 — 1)/2. For B^^^ = x < B* , 
the variables i?'-^-' decrease with r under iteration, and tend to zero for large r. For 



Linear polymers on Fractals 



9 





A® 



B 



(r) 



-(r) 



D' 



,(r) 



Figure 6. Restricted partition functions for the 3-simplex lattice. 




B 



B 
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Figure 7. Graphical representation of terms that contribute to the recursion equation for 
^(^+1) for the 3-simplex lattice. 



^(0) > 2?*, B^^^ increases with r, and diverges to infinity, and the series for P{x) and 
C{x) diverge. This implies that 

= B* = = [v^ - l]/2. (11) 

which shows that the growth constant fi of SAWs on this lattice is the golden mean. 
Using Eq. (jH)), we see that P{x) satisfies the functional equation 

P(x) = ^ + ^Pix' + x'). (12) 

From this equation, we get P{x = B*) = For x near B*, we can linearize the 

recursion equation. We write x = B* — 6, and Q{6) = P{B*) — P{B* — 5) giving 



Q((5)^/i-25 + ^g((2 + /i-2)5). 

If we assume that the singular part of Q{5) varies as this gives us 

log 3 



a = 2 



0.7342. 



(13) 



(14) 



log(2 + /i-2) 

The recursion equations for A'^'^^^^ and C'^'^^^^ are linear in A^'^'^ and C^^\ with coefficients 
that depend on B^'^\ If we start with a value 5 -C 1, then there is a value tq, such that 
for r — tq <^ —1, we can assume S^^^ ~ i?*, and for r — Tq ^ 1, B^'^''> ~ 0. As 
^(r+i) _ ^* ^ (2 + _ B*), (15) 
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Figure 8. Plot of the difference between the exact value of logP„ and the asymptotic 
form showing log-periodic oscillations with n for the 3-simplex lattice. 



we have clearly 
log(2 + /i ^) 



For r < ro, in the recursion equations for A^^^ and C'^^'\ we can put B^"^^ = B*, giving 

(17) 



This implies that A^^^ and C^^^ vary as A!j_, where A+ is the larger eigenvalue of the 2x2 
matrix in Eq. (fT7j) . For r > tq, 

^(r) ^ ^(ro)^ and - 0, then, in Eq. ©, the leading 
contribution comes from the term r = Tq, giving 

C(x = a;,-5)^(^ro ^(l/5)-7, (18) 
with 

= ^ 1.3752. (19) 

^ log(2 + /i-2) ^ ^ 

The typical size of the polymer for a; = — 5 is 2^o, which varies as {l/6y, with 

u = - — /"^^ - 0.7986. (20) 
log(2 + /i-2) ^ ^ 
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It is quite straightforward to generate the exact series P{x) using Eq.© with symbohc 
manipulation programs hke Mathematica. Just keeping terms up to r = 7, we get a series 
for P{x) exact up to x^^^ ( the coefficient of this term has more than 70 digits!). As 
seen above, the singular part of P{x) varies as (1 — x^Y""" for x near its n^^ Taylor 
coefficient must vary as fi^n°'~^. However, in general, the solution of Eq. lfT^ allows an 
additive term which is periodic function of log 5 with period log(2+/i~^). In Fig. |H1 we have 
plotted the difference A log P„ between the exact value of log -P„ and the asymptotic value 
log(/i"n°~^). We clearly see the log-periodic oscillations. Such log-periodic oscillations 
are seen in many other systems showing discrete scale- invariance [ 20^ . The existence of 
these oscillations makes it very difficult to estimate critical exponents from extraplations 
of exact series expansions using only a few terms [ |^ . 

For the fugacity x > Xc, the linear polymer fills the available space with finite density. 
Since the logarithm of the single loop partition function for the r-th order triangle is now 
proportional to the number of sites in the triangle, we define the free energy per site in 
the dense phase as 

log 5^'') 

fU) = lim . (21) 

Then, from Eq.®, it follows that f{x) satisfies the equation 

fix) = ^fix' + x'). (22) 
The density of the polymer p{x) in the dense phase is given by 

p{x) = x^fix) = ^T^Mx' + x'). (23) 
ax 6[x + 1) 

Iterating this equation, we get 

oo o 

P(^) = n[(^^^^ + o)(^^''^ + (24) 



From Eq. ()22|) . it can be seen that as x tends to Xc from above, the polymer density 
decreases as {x — Xc)^~". Equivalently, for small densities p, the entropy per monomer 
/i(p) ~ yu(0) — Kp^^ , where K is some constant. 

3.2. Quenched averages 

We note that the 3-simplex lattice is not homogenous, and different sites are not equiva- 
lent. In the context of disordered systems, the calculation of P{x) and C{x) are examples 
of annealed averages. We now show how one can calculate "quenched averages" in this 
problem. 

We define Pn{S) as the number of polygons of perimeter n that pass through a given 
site 5* ( these are called rooted polygons, one could also study rooted open walks), and 
ask how does it vary with S. What is its average value, variance etcetra ? A quantity of 
particular interest is the average value of logPn(5'). A good understanding of these for 
the regular fractals is prerequisite for understanding the more complicated case of random 
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Figure 9. (a) labelling of different sites on the 3-simplex (b) definition of weights for 
rooted graphs 



quenched disorder. We mention briefly some very recent results about the behavior of 
fluctuations of Pn{S) with S [l22]. 

To describe rooted polygons, we first need to set up a coordinate system to describe 
different sites on the fractal. A simple and natural choice for this is shown in fig. IHl A 
point on the r-th order triangle is labelled by a string of (r — l) characters, e.g. 0122201 . . .. 
Each character takes one of three values 0, 1 or 2. The leftmost character specifies in which 
of the three sub-triangles the point lies (0, 1 and 2 for the top, left and right subtriangle 
respectively). The next character specifies placement in the (r — l)-order sub-triangle, 
and so on. The restricted partition functions for the rooted polygons are also defined inlHl 
Here B^\x, S) is the sum over walks on the r-th order triangle that go through the left 
and right corners of the triangle, and also visit the site characterized by string S inside 
the triangle. Other weights are defined similarly. 

A site in the (r + l)-th order will be characterized by a string of r characters. Hence, a 
site characterized by string S at the r-th stage will be characterized by one of the strings 
OS", 15* or 25". The recursion relations in this case are linear, and can be written in the 
matrix form 

Bi^^'\x,sS) \=mA Bi^\x,S) , (25) 
\B(r^'\x,sS) J \Bj,^){x,S) J 

where 

/ B^ \ f B B^ \ / B B^ \ 

Mo= B B^ \ ; Mi = \ B^ \ ; M2= \ B^ B \ , (26) 

\ B^ B J \ B'^ B J \ B'^ J 

and we have suppressed the superscript (r) on B in the matrices. The generating function 
for rooted polygons is given by 

P{x,S) = J2B(:\x,Sr)B(^^\ (27) 

r=l 
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Figure 10. A plot of the average value of logP„(S')/i^" as a function of n. The uppermost 
and the lowermost curves are theoretically derived upper and lower bounds to this number 
over different positions of the root. The dashed and dotted lines show the annealed and 
quenched average values respectively. 



where Sr is the string consisting of the last (r — 1) characters of the position of S. If we 
ignore the constraint that the walk has to pass through S, we get an upper bound on 
the number of such walks, B^\x,S) < B^^\x), where inequality between polynomials is 
understood to imply inequality for the coefficient of each power of x. This implies that 
for all sites S, 

P{x,S) <Y.B'^'-^ . (28) 

r=0 

If we write the right-hand-side as then J-'u{x) satisfies the equation 

J^u{x) = x^ + J^u{x'^ + x^). (29) 

This equation should be compared with Eq.()13p. and differs from it only in the absence 
of the factor 1/3. This functional equation has the fixed point at x* = l//i, and linear 
analysis near the fixed point shows that J-'u{x) diverges as — log(l — xfi) for x tending 
to l//i from below. This implies that the coefficient of in the Taylor expansion of J-'u 
varies as fi^/n for large n. Thus, 

Pn{S) < Kjj'^/n, for all n, and all S, (30) 

where K is some constant. Sumedha and Dhar [|22] also derive a lower bound Pn{S) > 
for all n and all S, where b = 21og/i/log(2 + /i"^) ~ 0.92717, and Ki is 
some constant. For a randomly chosen S, the problem reduces to a random product of 
noncommuting matrices Mq, Aii and Ai2- The expected value of the logarithm of a 
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Figure 11. Different restricted partition functions for the 4-simplex lattice. The shaded 
squares denote graphs of r-th order 4-simplex, of which only the external bonds are shown. 



product of r independent random matrices varies linearly with r for large r. This implies 
that for large n, the quenched average < logP„(5') > varies as nlogfi — (2 — ag)log?7.. 
The numerical estimate of aq is approximately 0.729 ± .001, which is just a bit less 
than the annealed value 0.7342. In figure EH we have plotted the numerical values of 
the upper and lower bounds to log(P„(S')/i~"), and the annealed and quenched averages 
log < P„(S')/i~" >, and < log(P„(S')/i~") >. It should be noted that the lower bound is 
the best possible, in the sense that for each value of n, there is a finite density of roots that 
saturate the bound. The upper bound is not always optimal. Also, whether the annealed 
and quenched averages are very close or less so appears to be an oscillatory function of 
logn. 

4. RECURSION EQUATIONS FOR OTHER FRACTALS 

We now briefly indicate how the real-space renormalization group technique of the 
previous section can be extended to other fractals. 

4.1. n-simplex fractals for n > 3 

The case of n-simplex fractal for higher n is a straight forward extension of the technique 
used for the 3-simplex. The case n = 4 was studied in [EI? and the treatment was 
extended to n = 5 and 6 in [ "23^ . 

For higher n, the number of restricted partition functions we have to define to generate 
a closed set of recursion equations is larger. Here the permutation symmetry between the 
corner vertices of the n-simplex graph is very useful in reducing the number of variables 
needed. For n = 4 and 5, we need two functions A^"^^ and B^"^^ to generate the closed 
loops generating function P(x). For the open walks generating function C{x), four more 
variables are needed. Their definition is shown in Fig. ^2 They will be denoted by 
CW^^M^^W and F^'^\ 

The starting values of these weights are 

A(^) = x; = v^; P(^) = D(^) = E^'^ = F^'^ = 0. (31) 

The recursion relations for these weights are written by constructing graphs by all possible 



Linear polymers on Fractals 



15 



ways. This leads to [ ITT] 

^(r+i) = ^2 ^ + 2A^ + AA^B + QA^B^, (32) 

^(r+i) = A^ + AA^B + 22B\ (33) 

and similar equations for the other variables (omitted here). The expression for P{x) in 
this case is 

oo 

P(x) = ^(4AW +3AM ). (34) 

r=l 

There is a similar but more complicated expression for the generating function for open 
walks. We note that recursions of A^'^^^^ and 5^^+^) depend only on A^"^"* and B'^'^\ For any 
given value of x, we can determine P{x) numerically by explicitly iterating the recursion 
equations. One finds 

Xc = 0.4373, fi = 2.2866, (35) 

and the corresponding non-trivial fixed point is {A*,B*) = (0.4264,0.04998). Lineariza- 
tion of Eqs. and (jH^ around this fixed point gives largest eigenvalue Ai = 2.7965. 
We can express the end-to-end distance exponent u in terms of Ai using arguments similar 
to the case of 3-simplex. We get 

•■ 0.6740, (36) 



logAi 

and specific heat exponent 

a = 2-uDi = 0.6519. (37) 

To find 7, we consider the recursion equations for (7*^*"+^^ and D^'^'^^\ which are linear 
in C^^^ and D^^\ and diagonalizing the 2x2 matrix gives the eigenvalue A_|_ = 4.2069, 
which corresponds to 7 = log[(A^)/4]/ log[Ai] = 1.4461. 

For the 5-simplex lattice, the closed loop generating function can be found in terms of 
variables A^''^ and B^'^\ defined as before. In this case, the recursion equations are more 
complicated [I221- We list these here to show how the complexity of the polynomials rises 
rapidly with increasing n: 

^(r+i) ^ 132{B''+AB^)+A'^{1+18B'^+96B^)+A^{3+12B+78B^)+A^{Q+30B)+QA^, (38) 
^(r+i) ^ igg^s ^ 220AB^ + SSA^B^ + A^{1 + 13B + 2205^) + 2A\ (39) 

We omit the other equations. Starting with {A^^\ B^^^) = (x, 0), one finds a non-trivial 
fixed point for x = Xc = 0.336017, which corresponds to fi = l/xc = 2.97603 and the 
fixed point {A*, B*) = (0.3265, 0.02791). Linearization of the recursion relations Eq. ()38|) 
and (j39|) around this fixed point leads to one eigenvalue Ai = 3.13199 greater than unity. 
Using this value of Ai one gets z/ = 0.6049 and a = 0.5955. 

The exponent 7 is found from the recursion relations for the variables corresponding to 
a single end point. The largest eigenvalue A+ of the matrix which characterizes the linear 
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transformation of the recursion relations of these functions is A+ = 5.24398 corresponding 
to 7 = 1.4875. 

For the 6-simplex lattice, even for the polygon generating function, we need three 
restricted partition functions A^'^\ B^"^^ and C^^\ corresponding the cases where the walk 
enters and exits the r-th order subgraph once, twice or three times respectively. For 
the open walks, we need six more variables. For details, consult [ 23]. In this case, the 
nontivial fixed point is found to be 

{A*, B*, C*) = (0.262352, 0.017588, 0.0007011). (40) 

This fixed point is reached if we start with x = Xc = 0.27166, which corresponds to 
the connectivity constant fi = 3.68107. The linearization of recursion relations about this 
fixed point yields only one eigenvalue, Ai = 3.4965, higher than unity giving u = 0.5537 
and a = 0.5686. 

By diagonalizing the matrix corresponding to the linear recursions for weights corre- 
sponding to one end point, one gets the largest eigenvalue A+ = 6.26709. This gives 
7 = 1.500094 for the 6-simplex lattice. 

4.2. n-simplex lattice in the limit of large n 

One can explicitly determine the critical exponents for a few more values of n using 
a computer to generate the recursion equations. However, this soon becomes very time- 
consuming. If n is very large, some simplifications occur, and one can determine the 
leading behavior of the critical exponents of the SAW problem on the n-simplex in this 
limit [E3- discuss this below. 

We start by noting that as n increases, the probability of occurrence of loops in a 
random walk on the graph decreases, and as the loops become rarer, the random walk 
without self-exclusion should approximate the properties of the SAW. In particular, we 
would expect that the growth constant for SAWs on the ?7,-simplex should be approxi- 
mately equal to — 1. Let us denote, as before, the restricted partition functions for 
the r-th order lattice corresponding to polygon generating function by A^'^\ B^'^\ C^^^ . . ., 
corresponding to configurations where the SAW enters (and exits) the graph once, twice, 
thrice . . .respectively. Then for the fixed point corresponding to the swollen state, we 
would have A* ~ 1/n. Since B^"^^ < A^^'' , it is at most of order 1/n^, and similarly, C* is 
at most of order Thus, B*, C*, approach zero faster than A* as n is increased. 

It is straight forward to write down the recursion equation for A^^~^^^ if we neglect the 
terms involving B, C, etc. There is only one configuration of walk corresponding to the 
term A'^ ( we drop the superscript of A^*"^ for convenience), but (n — 2) terms of type A^, 
and [n — 2){n — 3) terms of type A^ etc. This gives us the recursion equation 

^(r+i) ^ ^2 ^ _ 2)^3 + _ 2)(n - 3)A^ + + {n- 2)!A". (41) 

When A is of order 1/n, each of these terms is of order 1/n^, and the sum needs to be 
evaluated with some care. In [^^, it was noted that, with only a small error when A is 
near A*, this series can be rewritten as 

j^ir+i) ^ A"(n-2)![l + l.A-^ + (l/2!)A-2 + ... + (i/(n-2)!)A-"+2] (42) 
~ (n-2)! A" exp(l/A) (43) 
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One can then study the asymptotic behavior of this recursion equation for large n. Here 
we shall use a different approach. If we change the variable from A^"^^ to e^^^ using the 
relation 

AM = -exp(^). (44) 
n \ n 



and use the fact that, for 3> r 3> 1, 

r 2 

n(l-j7n)^exp(-^), (45) 

the series of Eq.(35) can be written as 

1 " f(''V 

A<-'>.(^)i;exp(^--). (46) 

Substituting x for r / y/n and replacing the summation over r for ^/n ^ e^^'^ ^ 1 by 
integration over x from — oo to oo (for e ^ 1, the lower limit can be extended to — oo 
with negligible error) we get 

^^^+^^^(^)exp(^). (47) 

From this equation, we see that the nontrivial fixed point is given by 

e^") = e* ^ y^logn. (48) 

Let us now look at the equation for B'^''^^\ Again keeping only the terms involving A's 
alone, we get 

^C'^+i) ~ ^ 2(n - A)A^ + 3(n - 4)(n - 5)^*^ + . . . (49) 

Using the same approximation as before, we get 

D(r+i) / 1 \ V- I ^'^ V27r log n 

~ — ) ) rexp(^^ ~ — — , (50) 

where we have used the approximate fixed point value of e* from Eq. ()48|) . The fact that 
^(^+1) decreases faster than justifies neglecting these terms in determining the asymp- 
totic behavior of the recursions. The value of the derivative of the linearized recursion 
equation for A at the nontrivial fixed point is 

d 



^ /^log^_ (51) 



This implies that the critical exponent v is given by 
2 log 2 



V 



logn 



\ ^ — i 1_ higher order terms 

logn 



(52) 
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We note that the correction term to Eq. ()41|) involving B^^^ in leading order are of the 
form 

^^(r+i) ^ 2A^B{n -2){n-3)+ AA^B{n - 2){n - 3){n - 4) + .... (53) 

Using the estimate Eq.(50) for B^'''\ we see that near the fixed point, the error in Eq.(41) 
^^(r+i) -g q£ order n^^^"^. This implies that the fractional error in the value of u using 
Eq.(41) is of order n^^^"^. Thus Eq.(41) is much more accurate than Eq.(52), where the 
error is of order l/(logn). 

To calculate 7, one has to consider configurations with one endpoint of walk inside the 
graph. Again, keeping only terms involving powers of A, the recursion relation for the 
weight of configurations with the walk entering the graph once ( call it X) is found to be 

~ X[l + {n- l)A* + {n- l){n - 2)A*^ + .... + (n - l)!^*^"-^)]. (54) 

Using the arguments as before, this can be evaluated using the fixed point value of A* , 
giving 



~ Xy'nlogn. (55) 
It follows that, for large n. 



7 



1 '^f + higher order terms 

logn 



(56) 



Note that the critical exponents do not take mean-field values, even when the fractal 
dimension of the lattice becomes greater than 4. This is clearly because of the special 
structure of the n-simplex lattice, where, even though the fractal dimension becomes 
greater than 4 for large n, the spectral dimension remains below 2, and probability of 
intersection of paths of random walks remains large. Also that the non-analytical de- 
pendence of the type '"^^j^^" in the critical exponents on the lattice cannot be obtained 
from e-expansion like power-series expansions in deviation of dimensionality from some 
reference value. 

4.3. Modified rectangular lattice 

This lattice is interesting as its fractal dimension, and the spectral dimension are both 
rational numbers. Also, one can get its graph by selectively deleting some bonds from the 
graph of d-dimensional hypercubical lattice. Since the (r + l)-th order graph is formed by 
taking only two smaller graphs, the recursion equations involve only a small number of 
configurations, and are easy to write down. But the number of variables needed is larger, 
as the symmetry of the graph is lower than that of the HB{b, d) family. 

The behavior of SAWs on the d = 2 lattice was studied in The recursion equation 
for a polygon is written in terms of five restricted generating functions (see Fig. IT^ by 
constructing graphs by all possible ways [ [T7] . This gives 

^('^+1) = B{1 + D), B^'+^^ = A^ + C^ C(^'+') = 2AC, 

jj{r+i) = + 2DE, ^('^+1) = D\ (57) 
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Figure 12. Diagrams representing the restricted partition functions for the various ways 
the polymer can cross the r*^ order rectangle. 



The starting values of these recursions are 

The polygon number generating function P{x) is given by [ ITTj 



4 oo (U{r)\2 

n^) = T + E[^]- (59) 

^ r=l ^ 

Numerical iteration of these equations gives Xc = 0.5914 and fi = 1.6909. The fixed point 
occurs at 

{A*, B*, C*, D*, E*) = (0.5000, 0.4201, 0.4124, 0.1902, 0.0362), (60) 

giving one eigenvalue Ai = 1.6839 greater than one. Using this eigenvalue one finds 
u = log 2/2 log Ai = 0.66503 and a = 0.6699. A similar analysis for the open-walk 
configurations gives the critical exponent 7 = 1.4403. 



5. SELF-AVOIDING WALKS ON FRACTALS WITH DIMENSION 2 - e 

In the preceding sections we studied the properties of SAWs on several different fractals. 
For the calculation of critical exponents of SAWs on still other fractals, see for example [ 
It would appear that for each fractal, one has to write down the polynomial recursion 
equations, and calculate the values of critical exponents using the technique outlined. 
There is no simple expresion for the critical exponents as a function of the geometrical 
properties of the fractal, ( an improved Flory formula?) that would allow one to predict 
these without doing the full calculation. Unfortunately, this level of understanding of the 
problem is still not achieved. 

The best way to understand such systematics is to study the variation of exponents 
across a family of fractals as some property of the fractal is changed. For the n-simplex 
family studied earlier, the fractal and spectral dimensions did not tend to the same value 
for large n. We now discuss variation of critical exponents for SAWs on the GMb fractals 
family as the parameter b is varied from 2 to infinity. For large b, both the fractal 
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and spectral dimensions of the fractal tend to 2 where the properties of walks are well- 
understood. Such a study makes contact with the formal e-expansion technique which has 
played an important role in development of renormalization theory of critical phenomena. 

The self-avoiding walks for the GM2 fractal have the same exponents as for the 3- 
simplex. For small values of b, it is straight forward to generate the explicit recursion 
equations on a computer, and determine the critical exponents. These were worked out 
by Elezovic et al for b=2 to 8, and these studies were extended to 6 = 9 by Bubanja 
et al [i27j. The general form of recursion equations for the function A^^\ B^^\C^'^'* and 
D^^'^ ( definition of these is same as in Fig. ^ are of the form 



Air) 

C(r) 



(61) 



where fb, Pu, Pu, P21 and P22 are polynomials of Bi'^\ whose exact form depends on b. 

From analysis of these equations, it was found that the critical exponent u for SAWs 
takes the values 0.7986, 0.7936, 0.7884, 0.7840, 0.7803, 0.7772, 0.7744 and 0.77218 as b is 
varied from 2 to 9. It thus seems to converge to the Euclidean value i>2d = 3/4. However, 
the value of the critical exponent 7 changes from 1.3752, 1.4407, 1.4832, 1.5171, 1.5467 and 
1.5738 as b is varied from 2 to 7, and seems to diverge away from the known exact two 
dimensional value 72^ = 43/32 ~ 1.344. 

The variation of these exponents with b for large b was explained in [|2H1 using the finite 
size scaling theory. For large b, the growth constant of SAWs on the 6-fractal would be 
close to the critical value in two dimensions. It is convenient to change variables from 
B^^^ to a variable that is proportional to the departure from criticality in these systems. 
We write 

5M=5,*_exp(eM). (62) 

Then is related to the properties of SAWs that traverse an equilateral triangle of 

side b from one corner to another. From finite-size scaling theory, this would be expected 
to be a function of a single variable eb^^'^: 

/,(i?)^^exp[^(66i/^)], (63) 

where K is some constant, and a is an exponent. From the conformal field theory [IHll29j. 
the scaling dimension of a spin at the corner of a wedge of angle 7r/3 is known, and that 
implies that a=15/4. The scaling function g{x) has to have the following properties: 
(i) It is a monotonically increasing convex function of x. We can set g{0) = 0, by redefining 
K. (ii) For e < 0, /& should decrease exponentially with b. This implies that g{x) ~ — la^l*^ 
for X <^ —1. (iii) For fixed e > 0, /b should vary as exp(6^) for large b. Thus, we must 
have 

g{x) ^ Kix^", for x > 1. (64) 

From Eqs. ()61|) . it is easy to see that the fixed point value of e for large b is given by 
g{elb^/'') ^ a log 6. Using Eq. (jM|) . this gives 

(66) 
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The derivative of with respect to B^^^ at the fixed point is found to be 

Ai(6) = b'/''^U=et ^ K^b'/^ilogbf-^, (66) 

where K2 = 2iJa{Ki/ a)^/'^^ is a constant. Expressing the critical exponent v{h) in terms 
of Ai(6), we get 

1 = 1 , 2z/- lloglogfe 
v{h) V 2v log6 ' ^ ' 

where the dots represent terms of order 

It is interesting to note that the leading correction term to the asymptotic value of v 
is negative. Thus this analysis predicts that variation of v with h is not monotonic. In 
particular, z/^ should be less than v for large enough h. This rather striking prediction 
of the finite-size scaling analysis given above was checked by Milosevic and Zivic using 
numerical Monte-carlo renormalization group studies [ EHl 1^ ■, and they found that this 
happens for h ~ 27. 

We can similarly determine the leading correction to the susceptibility exponent 7. 
It is known that in two dimensional critical theory, a 3-leg vertex has a higher scaling 
dimension than a 1-leg vertex. This implies that at the critical point e = 0, the ratio 
C(^)/AW will tend to zero as a power of 6, and hence C^"^' can be ignored in determining 
the leading correction. Then the value of the exponent 7 is determined by the value of 
P\\ at the fixed point e = e*(&). We can now write the scaling ansatz for this variable: 

Pn(e,&)^A^3&^exp[/i(e6i/-)], (68) 

where K^^ is a constant, c is a critical exponent, and h{x) is a scaling function. 

For e < 0, it is known that as 6 ^ cxd, varies as [12^1 E21- This implies that 

we must have h{x) ~ (—1/64) logx, and c = 1/48. 

Now, for large x, the functions h[x) and g[x) should increase in a similar way. Then 
using Eq. ()63|) . we can write P\\ as 

Pn(eb,&) ^ K,W^''^^^\h{tlh^l-)-g{tlh^l-)\B\h). (69) 

For large a;, both h[x) and g[x) increase as K\x^^ ^ but the leading dependence is the same 
exactly. Thus, for large x, h{x) — g{x) varies at most as a sublinear power of log 6. Thus, 
we get 

lim — ; = c -I- a. (70) 

fe^oo log 

As 7(6) = log[ ^^^^^^ ]/ logfe, this implies that 

lim 7(6) = 2(c + a- l)z/ = 133/32. (71) 

Using the known scaling exponents for the dense polymer problem in two dimensions, it 
was shown in [ |2H] that the leading correction to the large h limiting value of 7 is given 
by 

7(6) = "^5 "^f ^ + higher order terms. (72) 

' 32 128 logfo ^ ^ ' 
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We note that the leading correction to asymptotic value of the exponent is proportional 
to 2 — df,. The next correction term is of order 1/ log 6, and is proportional to 2 — D^, where 
Df, = 2 — e is fractal dimension. These are like the e-expansions, except that there are 
several inequivalent definitions of dimension for fractals. Thus there are several e's, and 
the exponents on fractals may require a multi- variable e-expansion. Interestingly, at higher 
orders, corrections to scaling to the finite-size scaling functions f{x),g{x) etc. would give 
corrections to the exponents of the type These are of the type exp(— A/e), and 

such correction terms are not calculable within the conventional e-expansions framework. 

The reason the critical exponents in the large b limit do not tend to the Euclidean 
value may be understood as a crossover effect. For large b, the space looks Euclidean at 
length scales smaller than b, and the effective polymer exponents for e > 1/b^^'^ would be 
near the Euclidean values. However, for e -C 1/b^^'^, the polymer has to go through the 
constrictions, and the asymptotic value of exponents for large polymers can, and do, take 
different values. 

6. THE COLLAPSE TRANSITION IN POLYMERS WITH SELF- ATTRACTION 

We have seen that the qualitative features of the behavior of linear polymers with 
no interaction other than the excluded volume interaction is well-modelled by SAWs 
on fractals. Now we will like to show that polymers on fractals can also be used to 
understand more complicated behavior of polymers like the transition from the swollen 
state to compact globular state transition that is observed in dilute polymer solutions as 
the temperature is lowered [1^. 

In order to model this, we have to include the attractive interaction between different 
parts of the polymer when they are close by, but not overlapping. The simplest lattice 
model for this phenomena associates an energy —E^ for each pair of nearest-neighbor 
lattice points occupied by the polymer that are not consecutive along the walk. The 
equilibrium properties of such a self-attracting SAW can be described using the grand 
partition function G{x,u); 

G{x,u) = J2 ^{N,Nu)x^u^-, (73) 

N,Nn 

where Q{N, N^) is the number of different configurations per site of a self-avoiding polygon 
of N steps and energy —NuEu, u = exp{PEu), x is the fugacity per step of the chain, and 
f3 is the inverse temperature. 

For small (3, the effect of Eu can be ignored, and the typical size of polymer varies as 
A^"^, where the exponent u takes the value for SAWs. This is called the swollen state of 
the polymer. For large (3, the polymer folds up like a tangled ball of yarn, in order to 
minimize the energy —N^Ey^. In this phase, called the collapsed phase, the typical size of 
polymer varies as N^^*^ in d dimensions. 

In the limit when the number N of monomer units goes to infinity, there is sharp 
transition from the swollen to collapsed phase at a critical value of u. This transition is 
described as a critical phenomena analogous to a tricritical point for magnetic system [|2]. 
For large polymers, the average gyration radius R at the transition behaves as ^ N"'^ 
where the exponent uq is intermediate between the value u for swollen state and the value 
i>c = ^/E> for compact globule on a lattice of fractal dimension D. 
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6.1. Self-interacting polymer on the 3-simplex lattice 

The properties of a polymer chain with self-attraction on a fractal were first studied by 
Klein and Seitz [E^- They used the self- avoiding walks on the Sierpinski gasket, which 
is the 6 = 2 member of the Given- Mandelbrot family. We consider below the case of 
3-simplex, which is somewhat simpler to treat. 




Figure 13. Restricted partition functions for the self-attracting walk for the 3-simplex. 



To calculate G{x, u) for the 3-simplex lattice, we define restricted partition functions 
Bq^ and Bi^\ for walks that cross the r-th order triangle, as shown in Fig. ^1 Bq^ is 
the sum of weight of all walks that enter an r-th order triangle of the 3-simplex from one 
corner and leave from another corner, but do not visit the third corner. B^^ is the sum of 
weights of walks that enter and leave the r-order triangle, and also visit the third corner 
of the triangle. Then it is easy to see that these weights satisfy the recursion equations 

= [Bo + B-,f + Bo[Bl + 2So5i + Blu], (74) 

= B^[Bl + 25o5i + Blu]. (75) 

The generating function for all loops is given by the formula 

oo 

PM = Y.[Bt^ + Bt^f?>-~\ (76) 

r=l 

We start with the initial weights 

B'i^ = 0; 5}'^ = X. (77) 

These variables tend to the trivial fixed point Bq = Bl = 0, if the starting value of x is 
less than a critical value Xc{u), and to the fixed point Bq = Bl = oo, if x > Xc{u). For 
X = Xc{u), b'(^ tends to zero for large r, and the Eq. (f?^ reduces to Eq.®. The critical 
properties of the chain are continuous functions of u, and there is no phase transition as 
a function of u. 

We note that here the recursion equations involve the interaction parameter u. This 
complicates the analysis. Consider a modified problem where the interaction —U occurs 
only between the nearest neighbor bonds in the same 2-nd order triangle, and not oth- 
erwise. One would expect the qualitative behavior of the system very similar, but now, 
the recursion equations do not involve u. In fact, they are the same as the case with- 
out self-attraction [Eq.®]. This observation is helpful in studying the collapse on other 
fractals. 
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Figure 14. Different fixed points for tlie 4-simplex lattice. Tlie non-trivial fixed points 
labelled C,T and 5* correspond to the collapsed, the tricritical 6'-point, and the swollen 
state. 



6.2. Collapse transition on n > 3 -simplex lattices 

The collapse transition on the 4-simplex lattice was first studied in [EH]- We restrict 
the attractive interaction to nearest neighbors within the same 2nd order simplex. Then, 
the recursion equations Eq. ()32H33|) describe the collapse transition also, if we use the 
starting weights 

^(2) =x^ + 2x^u + 2x\^, and S^^^ = x\^. (78) 

The renormalization group flows of the variables A and B in this case is shown in 
Fig. EI There is the trivial fixed point A* = B* = 0, which is reached if x < Xc{u). 
If X > Xc{u), the fixed point A* = B* = (yo is reached. The two-dimensional space of 
possible initial conditions (A^^\ B^^^) is divided into the basins of attraction of these two 
fixed points. The common boundary of these basins is a line. This line is an invariant 
sub-manifold of the renormalization flows (i.e. points starting on the line remain on the 
line). On this line we have three fixed points : 

1. The fixed point (0.4264, .04998) corresponds to the swollen state. This is reached 
for u < Uc — 3.106074 and x = Xc{u). This is denoted by S in Fig. El If we start 
at a point near this fixed point on the critical line, the renormalization flows are 
towards this fixed point. 

2. The fixed point (A*, B*) = (0, 22"^/^) is reached when u> Uc- It gives one relevant 
eigenvalues A = 4 corresponding to z^c = I/-D4 = 1/2. In this phase, the polymer fills 
the available space with a finite density of monomers. This fixed point is denoted 
by C in Fig. El This is also attractive for points starting nearby on the critical 
line. 

3. The fixed point (1/3,1/3), denoted by T in Fig. HH is obtained for u = Uc- On 
the critical line x = Xc{u), this unstable fixed point lies between the fixed points 
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corresponding to the swollen and the collapsed phases. The matrix corresponding to 
linearized renormalization transformation at this point has two eigenvalues greater 
than one; Ai = 3.7037 and A2 = 2.2222. The third relevant eigenvalue corresponds 
to the renormalization of fugacity of end points ( A+ in Eq. (18)). Thus this is a 
tricritical point, at this point, = ln2/lnAi = 0.5239. The singular part of the 
free energy varies as {u — Uc)^~° with 

a = 2- = 0.36027. (79) 
In A2 

One can similarly study the n = 5 simplex lattice [|^ , with the interactions are confined 
to the internal bonds of the second order graph. It is found that just like the n = 3 case, 
there is no collapsed phase. One can understand this by noting that the ground state 
configuration corresponds to a walk that visits all the 5 sites of the 2nd order subgraph. 
There are many configurations corresponding to the minimum energy, and most of these 
are extended. Even when there is attractive interaction between all pairs of neighbors, 
the energy-cost of pulling a polymer confined to an r-th order subgraph to something that 
is spread over a subgraph of order (r + 1) is finite. As there is a large entropy associated 
with the place where this break occurs, this makes the collapsed phase unstable to such 
breaks, and brings the collapse transition temperature to zero. 

In a similar way, we can study the 6-simplex. Here we have coupled polynomial re- 
cursion equations of degree 6 for the three variables A^'^\ B^"^^ and C^^\ For the details 
of the recursion equations, see [ |21] . Again, we have two trivial attractive fixed points 
corresponding to {A*, B*, C*) equal to (0, 0, 0) or {00, 00, 00). There is a two-dimensional 
critical manifold that marks the boundary of the basins of attraction of these fixed points. 
For the renormalization flows on this two-dimensional critical surface, there are two at- 
tractive fixed points: 

The fixed point {A*, B*, C*) = (0.262352, 0.017586, 0.000701) corresponds to the swollen 
state as discussed in Sec. 4.1. 

The fixed point {A% B% C*) = (0, 0, 0.071329) is reached for all u > Uci= 3.4999847) at 
X = Xc- The largest eigenvalue for this case is A = 6 corresponding Uc = 1/ Dq = 0.3869. 
This describes the collapsed phase of the polymers. 

The common boundary of basins of attraction of these two fixed points is a 1- dimen- 
sional line, which is also an invariant submanifold for the renormalization flows. This line 
has one attractive fixed point: 

The fixed point (0.12949,0.09572,0.05344) is obtained for u = and x = Xc and has 
two eigenvalues greater than one; Ai = 5.4492 and A2 = 1.9049. This is a tricritical point 
with exponent z/j = 0.4088 and a = 2 - In A1/A2 = 0.6309. 

The crossover exponent at the tricritical point is = Vtjvth = 0.3801, where the expo- 
nent Uth controls the divergence of thermal correlation and is defined as Uth = In 2/ In A2 = 
1.0755. 

There are two other fixed points (0.254037, 0.022159, 0.07098) and (0.2000, 0.0666, 
0.0666) on the line separating the basins of attraction of the fixed points corresponding to 
the collapsed and the swollen phases. These are purely repulsive, and cannot be reached 
starting with our choice of initial condition. These correspond to higher order multicritical 
points( tetra-critical). 
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We note that the collapse transition on fractal lattices corresponds to a new fixed point, 
intermediate between swollen (SAW) phase and the collapsed phase and cannot be viewed 
as a perturbation of the Gaussian fixed point describing random walks. 

6.3. Some unusual phases 

The nature of possible collapsed phases of linear polymers on fractals with a finite 
ramification number depends strongly on the geometrical constraints imposed by the 
connectivity properties of the the fractal, much more so than in the extended phase. 

As we already noted, for the n-simplex lattice with n odd, a collapsed phase with a finite 
density of the polymer is not found. For the modified rectangular lattice, the behavior 
of linear polymer with self-attraction was studied in [ [33] . Here the number of variables 
needed to describe the polymer with self-attraction becomes rather large. For example, to 
describe closed loops it is necessary to introduce additional weig 

where the subscript 1(2) indicates the number of extra corner sites of the r*'* -order are 
visited, a total of eleven variables. To describe open chains, we would need 17 additional 
variables, making a total of 28 variables- a rather formidable number. However, as in case 
of ra-simplex lattices many of these variables are irrelevant and may be set equal to zero. 
Equivalently, we study collapse when the attractive interaction is restricted to the first 
order rectangles. If we restrict ourselves to the analysis of closed polygons, we require 
only the five variables defined earlier in section 4.3. 

Interestingly, one finds three different phases of the polymer, depending on the value of 
the attraction parameter u. For small u < Ud = 3.2023, the polymer exists in the swollen 
phase, with typical size R ~ N'^ , and u = 0.6650. For large u > Uc2 = 3.2341, it exists 
in a compact phase of finite density, with R ~ iV^/^. However, between the swollen and 
compact phases, for Ud < u < Uc2, one finds a 'rod-like phase', where the shape of the 
polymer is highly anisotropic. In the x-direction, the average extent of polymer increases 
as A^, and in the other direction it remains finite. At u = Ud, there is a non-trivial 
fixed point of the recursion equations of period 2. Linear analysis of the renormalization 
equation about this fixed point gives u = 0.80503. At u = Uc2, one gets u = 1/2, same 
as in the collapsed phase, with possible logarithmic corrections. There is a logarithmic 
singularity in the specific heat also. The reason why the rod-like phase is found in this 
case is not clear : obviously the anisotropy of the lattice is responsible for it, but it does 
not matter so much in other( swollen or collapsed ) phases, where the anisotropy is of 
the usual type ( the ratio of average diameter of an A^-stepped walks in the x- and y- 
directions is a finite number). 

Another different type of phase, labelled 'semi-compact' phase was found by Knezevic 
and Vannimenus in their analysis of the collapse transition on the HB{3, 3) fractal [I37j. 
In this case, the connectivity of the graph is such that a linear polymer cannot fill the 
available space with a finite density. For large value of the attraction strength u the 
polymer shrinks into a "semicompact" state. In this phase, the average monomer density 
tends to zero for large polymers. 

The recursion relations found for restricted generating functions A^'^^ and B^^^ where 
A^"^^ counts the number of configurations when the polymer goes once through the r*'^- 
order gasket while B^^^ counts the configurations where the polymer goes twice through 
the gasket, are rather complicated [EI], but for the large u regime, the equations are 
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dominated by only a single term in each polynomial. For large order of iteration, the 
recursion equations in this phase may be approximated by 

^(r-+i) ^ 320^^5^; B^''+^^ ~ 4308^^5^ (80) 

These can be linearized by taking the logarithm of both sides, which shows that for large 
r, logA^^^ and log 5^^^ vary as A!^, where X± = (11 ± v^73)/2. They will both diverge 
to +00, or to — oo depending upon the coefficient of proportionality being positive or 
negative ( both are of same sign for the largest eigenvalue). Thus, this constant must be 
proportional to distance from the critical line x — Xc{u). If Xc — x = 6 <^ Xc, number of 
iterations before the deviation becomes of order 1 is 

ro^log(l/5)/logA+. (81) 

As the size of the polymer is approximately 2'"° ~ (1/5)*^', this corresponds to u' = 0.48195. 
The fractal dimension of the chain in this phase is, l/u' = 2.07491 a value just a bit less 
than the fractal dimension of the lattice D = 2.09590. For x = Xc, A^^^ decreases to zero 
as exp(— aA^) where a is some constant, and B^'^^ increases to infinity as an inverse power 
of A^'\ 




Figure 15. Schematic representation of the different phases of a linear chain with self- 
attraction on the Sierpinski gasket with at most two visits per site allowed: (a) the swollen 
phase (b) branched polymer phase with branches made of the doubled-up chain (c) the 
collapsed phase 



In all the cases discussed so far a restriction on the walks that a lattice bond or a lattice 
site visited once cannot be visited again has been imposed. One can relax this condition 
and allow the walk to revisit a site already visited so long as no bond is traversed twice. 
A consistent set of visited bonds is called a trail, if different sequences in which the 
same set of bonds bonds may be visited are considered equivalent [ |SH] • Orlandini et 
al [ESI considered self-attracting SAWs on a 2d Sierpinski gasket, with this constraint, 
and showed that this model displays a new multicritical point corresponding to a collapse 
from linear into branched polymer in which the polymer becomes a randomly doubled-up 
chain ( see fig. [TS|) . which is followed by a further transition into compact globule. In 
Fig. ^[ we have drawn schematically what the different phases look like ( drawn here for 
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the Euclidean 2- dimensional space). The universality class of the two new multicritical 
points corresponding to the phase transition points is different from the tricritical points 
discussed so far. 

7. BRANCHED POLYMERS 

The models of branched polymers is related to other important problems in statistical 
mechanics, such as lattice animals, percolation [HO] and the Lee- Yang edge singularity [ 
1^ . The study of branched polymers on fractal lattices are therefore very instructive as 
one can analyze the system in detail. The number of different branched polymers made of 
N monomers and their average gyration radius are expected to grow as fj,'^N~^ and N", 
respectively for large A^; 6 and u are critical exponents, whose values are different from 
those from that of linear polymers. If we allow loops in the cluster, we get the model of 
lattice animals, which corresponds to the p —>■ limit of the percolation problem. These 
are known to be in the same universality class as branched polymers. If in addition to the 
excluded volume interaction, one has an attractive nearest-neighbor interaction between 
the monomers, the branched polymers can undergo a collapse transition, just like the 
linear polymers. Near the collapse point, the free energy per monomer shows singular 
critical behavior. 

Knezevic and Vannimenus realized that real-space renormalization method for studying 
linear polymers on fractals can be extended directly to the case of branched polymers [ 
1321 IS] • They considered asymptotic properties of branched polymers with attractive self- 
interaction on fractal lattices, restricting the attractive interactions to bonds within first 
order units of the fractal lattices. We summarize their findings here. 

7.1. The 3-simplex lattice 

The macroscopic thermodynamic quantities of interest can be obtained from the grand- 
partition function G{x, u), whose definition is same as Eq. (f7H|) . with Q{N, Ny) now defined 
as the average number of unrooted branched polymers with bonds, and A^^ nearest 
neighbor bonds, the average being taken over different positions of the polymer on the 
fractal. This can be determined in terms of the r-th order restricted partition functions 
defined in Fig. UHl 

The closed set of recursion equations involve six restricted generating functions as 
described in Fig are easily written down [ |32] : 
^{r+i) ^ A[1 + 2B + 2B^]+2B^C + F[B^ + + 2BD], 

^(r+i) ^ B"^ + B'^ + F[ABC + 2AB] + F'^[B + D], (82) 

and similar equations for other variables. G{x, u) can be seen as a sum of terms of the 
same general form as Eq. (0) and (jH)). The singular behavior of the sum can be analyzed 
by looking at the fixed points of the recursion equations. There are two trivial fixed points 
with all variables zero or infinite. For any given value of m, we have to tune the initial 
value of X to a critical value Xc{u) to get a non-trivial fixed point. 

The analysis of the recursion equations lead to following three different nontrivial fixed 
points: 

1. For all u < Uc = 5.5, and x = Xc{u), the recursion equations lead to a non-trivial 
fixed point {A*, B*, . . .) corresponding to the swollen state with large scale properties 
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Figure 16. Diagrams representing the six restricted generating functions for branched 
polymers on the two dimensional Sierpinski gasket. C^*"-* corresponds, for instance, to 
configurations where a part of the polymer joins two vertices of an r*^ order triangle while 
one of its ends penetrates through the third vertex. The diagram on the right shows a 
term B^'-^C^''^ F^'-'> contributes to B'^'-^^l 

same as for the random-animal problem (m=1). All the fixed point values A*, B*, . . . 
at this fixed point are non-zero. Linear analysis of the recursion equations about this 
fixed point shows that the generating function G{x, u) has a power-law singularity 
when X = Xc{u) and the critical exponents are u = 0.71655 and 6 = 0.5328. 

2. For u > Uc, and x = Xc{u), the polymer is in the collapsed phase. The corresponding 
fixed point has A* = B* = D* = 0. F^"^) tends to zero, and C^'^^F^'') and E^"^ F^"'^'^ 
tend to a finite limiting value. The largest value of the linearized renormalization 
transformation matrix is 3, corresponding io Uf. = l/D^ = 0.63093. This dense 
phase of branched polymers is in the same universality class as the spanning trees. 
For the spanning trees, one can also define a chemical distance exponent z, by the 
relation £ R^, where i is the distance along bonds between two randomly picked 
points on the polymer, and R is the euclidean distance. This was calculated in [ 
US], and shown that z = log[(20 + y205)/15]/ log2 ~ 1.1939. 

3. The collapse transition occurs at m = Uc- At this fixed point, all the values A*,B*, . . . 
are non-zero. This is a tricritical fixed point, with two eigenvalues larger than 
one. The exponent Ut = 0.63250 is very close to z/^. The exponent a is negative, 
a = —4.0269, showing that the singularity in the specific heat at u = Uc is very 
weak. 

The closeness of Ut to z/^ has also been found on square lattice where an accurate transfer 
- matrix study [|33j of the collapse of branched polymer gives z/j = 0.509 ± 0.003 which is 
very close to the compact value 1/2. This suggests that the phenomenon is not accidental 
and should have a more general explanation. 

The problem of linear polymers is recovered if all terms containing the functions E and 
F one suppressed in Eqs.(7.1) and (7.2). The truncated equations have only one fixed 
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point, corresponding to the swollen phase with u = 0.7986 and there is no collapsed phase 
(see Section 3). 

7.2. The 4-simplex lattice 

The closed recursion equations in this case involve eleven restricted generating functions. 
The number of polymer configurations to consider is large, and Knezevic and Vannimenus 
[131] used computer enumeration to sort them out. We omit the details. 

The analysis found three non-trivial fixed points that describe the large-scale behavior 
of self-interacting branched polymers. 

1. For M < Mc ~ 2 the random - animal fixed point corresponding to the swollen phase 
of the polymer is reached. Linearizing around the fixed point one finds only one 
relevant eigenvalues Ai = 3.14069. With this eigenvalue one finds u = log 2/ log Ai = 
0.60566, and 9 = 0.75667. 

2. For u > U(. one gets a fixed point corresponding to the collapsed state in which 
polymer occupies all vertices of the lattice. The relevant eigenvalue found in this 
case is A = 4 which gives Vc = 1/2. 

3. For u = Uc the fixed point reached represents a tricritical point as it has two 
eigenvalues greater than 1: Ai = 3.94050 and A2 = 1.32094. The exponent Uf = 
log2/logAi = 0.50546. This value is very close to the value Uc = 1/2 of the col- 
lapsed phase. The exponent a = —2.8267 is negative but less negative than for the 
3-simplex lattice. 

7.3. Other fractals 

Knezevic and Vannimenus also studied the branched polymer problem on the GM^ 
fractal [ |32l IHI • In this case, they found the unexpected result that unlike the case on 
the GM2 fractal ( for which the behavior is the same as on the n = 3 simplex), for the 
6 = 3 case, the number of animals of size n grows as /i" exp(Kri'^), where < ip < 1. This 
corresponds to an essential singularity in the generating function of branched polymers 
G{x,u = 1) ~ exp(|^-^), with p = tfj/il - tjj). 

The analysis of Knezevic and Vannimenus has recently been extended to all b, and one 
finds that for b > 3, the average number of animals per site behaves as /i(6)" exp(fi;n'^), 
where the values of the the singularity exponent ip, and the size exponent u can be 
determined exactly [EE]. 

8. SURFACE ADSORPTION 

It is well known that a long flexible polymer in a good solvent can form a self-similar 
adsorbed layer near an attractive wall at the critical temperature T^. Using the correspon- 
dence between an adsorbed polymer chain and the model of ferromagnets with n-vector 
spins in the limit n ^ with a free surface, it has been shown that the adsorption point 
Ta corresponds to a tricritical point and in its proximity a crossover regime is observed. 
In particular, the mean number of monomers, M, at the surface is shown to behave as [ 

im 
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M = {Ta- tY^/'^-^\ for T < Ta] 

= iV^ forT = T,; 

= (T-T„)-\ for T>T,. (83) 

A good model for this phenomenon is a SAW on some lattice with an absorbing surface 
(boundary). Every site on surface visited by the polymer contributes an energy —Eg. This 
model has widely been studied on various lattices and via a number of techniques that 
include exact enumeration [11111321; Monte Carlo [EOj; transfer matrix [Ej, renormaliza- 
tion group [|221 etc. For a 2-d Euclidean lattice, exact value of found from conformal 
field theory is 1/2 [E^. 

Bouchaud and Vannimenus [EH were the first to apply RSRG techniques on fractals 
to study the linear polymers near an attractive substrate. They showed that the known 
phenomenology of the adsorption -desorption transition is well-reproduced on fractals, 
and different critical exponents can be evaluated exactly. The values of the exponent (p 
for HB{2, 2) and HB{2, 3) fractals were found to be 0.5915 and 0.7481 respectively. They 
also showed that for a container of fractal dimension D and adsorbing surface dg, has 
lower and upper bounds; 

l-{D-ds)u<<f)<^. (84) 

8.1. Surface adsorption in polymers with self-attraction 

Polymer chains with self attraction near an attractive surface can exhibit a rich va- 
riety of phases, characterized by many different universality domain of critical behavior 
[57]. This is due to competition between the two interactions. At the intersection of two 
tricritical surfaces, one corresponding to the 6'— transition and another to the adsorption- 
desorption transion, one can expect higher order critical points. Understanding of the dif- 
ferent phases that are possible and understanding and classifying the multicritical points 
that can exist appear difficult on standard Euclidean lattices[47]. It is here fractal lattices 
have been particularly helpful [58]. The values of critical exponents found are, of course, 
different for different fractals, but the general features of the phase-diagrams remain the 
same. Investigating the problem on fractals helps us understanding the problem in real 
experimental systems. 

We consider a linear polymer chain on a truncated n-simplex lattice and make one 
boundary surface of it attractive. We treat one of the edges of the fractal container as 
an attractive surface and associate an energy —Eg < with each site on it occupied by 
the polymer, and an energy Et > 0, with each occupied site that is at a distance 1 from 
the occupied surface i.e. on the layer adjacent to the surface, and an energy —E^ for 
nearest-neighbor bond between monomers not consecutive along the chain. 

Then, to each iV-step walk having A^^ steps along the surface, Nt steps lying on the 
layer adjacent to the surface and A^^ number of nearest neighbors a weight u^^t^^u^^ 
is assigned, where uj = exp{PEs), and t = exp{—PEt). The grand partition function for 
this system is given by 

G{x,uj,t,u)= J2 ^{N,Ns,Nt,Nu)x^uj^H^'u^\ (85) 

N,Ns,Nt,Nu 
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where Q{N, Ns, Nt, Nu) is the number of different configurations per site of a S'AVT of 
steps, rooted at a specified site on the attractive surface, with given values of Ng, Nt and 
N^. We may put Et = (or t = 1) for simphcity, but a non-zero value seems to give rise 
to interesting behavior, and is also important to many physically realizable cases. The 
attractive interaction between monomers is restricted, as in preceding sections, to bonds 
within the lowest order subgraphs of the fractal lattices. 

As shown in section 6.2, polymers in the n-simplex container show a collapsed phase 
only for even values of n. Thus, the 4- and 5-simplex lattices exhibit contrasting behavior 
and represent two different scenario which may arise in real systems. In the case of 5- 
simplex lattice, there is no collapse transition possible in the bulk, but its 4-simplex surface 
can show a collapsed phase. In the case of 4-simplex lattice, on the other hand, there 
is a collapsed globule phase in the bulk but no collapse possible in the surface- adsorbed 
polymer. We can therefore have situations in which the bulk acts as a poor-solvent 
medium, while the surface acts as a good solvent medium or the opposite case of surface 
being poor-solvent medium and the bulk good solvent. However, the situation in which 
both bulk and surface show collapsed phases can not be modelled by polymers on an 
77,-simplex fractal. 

8.2. The 4-simplex lattice 




Figure 17. The restricted partition functions for the rth order 4-simplex with attractive 
interactions at one boundary surface. The internal structure of the 4-simplex is not shown. 
The shaded triangle reprsents the attractive surface. 



For the 4-simplex lattice the grand partition function of Eq. ()85p can be written in 
terms of the five restricted partition functions shown in Fig 17. The shaded regions in 
the figure represent the surface. Out of five configurations, two (A^"^^ and B^^^) represent 
the sum of weights of configurations of the polymer chain within one r-th order subgraph 
away from the surface, and the remaining three {S^'^\ C^'^^ , and E^'^^) represent the surface 
functions. The recursion relations for these restricted partition functions can easily be 
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Figure 18. Sections of the oo — u — t phase diagram for two different values of t for the 
4-simplex lattice: (a) t= 0.2 (b) t= 0.5. Regions marked by AS, DS and DC represent, 
respectively, the adsorbed polymer in swollen state, desorbed polymer in swollen and col- 
lapsed (globular) state. The collapse transition between the DC and DS phases is denoted 
by the dashed line. The special adsorption line is indicated by full line and part of it by 
dotted line. The point where 6'-line meets with the adsorption line is a multicritical point. 
The dotted part of the adsorption line indicates the region of coexistence of adsorbed 
SAW and the (bulk) globule phase. 



written (see [58] for details). The equations for A^^^^'' and 5^''+^) do not get affected by 
the presence of the surface interactions. 

In this case, we have three phases possible: the desorbed swollen(DS), the desorbed 
collapsed (DC) and the adsorbed swollen (AS) phases. The fixed points corresponding 
to the desorbed phases have S* = C* = E* = 0, and A*,B* equal to the value for the 
4-simplex at the S or C fixed point ( section 6.2). The fixed points corresponding to the 
polymer adsorbed on the surface, on the other hand, have A* = B* = C* = E* = 
and S* equal to the value corresponding to the swollen phase fixed point on the 3- 
simplex surface ( section 3.1). The phase boundaries, determined numerically by finding 
the basins of attraction of these fixed points are shown in Fig. ^| 

If we start at any point on the boundary between the DS and AS phases, under renor- 
malization, we flow to the symmetrical fixed point S* = C* = A* , B* = E* , with A*,B* 
having the value corresponding to the point S in Fig. Linearizing the recursion 

equations near this fixed point two eigenvalues A = 2.7965 (corresponding to the swollen 
bulk state) and A,^ = 2.1583 greater than one [58] are found. This point is the expected 
symmetrical special absorption point which describes the polymer at the desorption tran- 
sition. A simple calculation gives the crossover exponent = 0.7481 and a = 0.6653. In 
Fig 18, the full line shows intersection of this surface with the surface t = constant. 

The bulk transition between DC and DS phases is described by the fixed point S* = 
C* = E* = 0, and A* and B* have the value (1/3,1/3) corresponding to the point T in 
Fig. 1131 This boundary has the equation u = Uc, and the value Uc does not depend on u 
or t. 
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The points on the boundary between the DC and the AS phases are found to fall in the 
basin of attraction of two different fixed points: FPi = {A* = S* = C* = 0, B* = E* = 
0.3568) and FP2= {A* = C* = E* = 0,B* = 0.3568, S* = 0.61803). Correspondingly, we 
see two different behaviors of various quantities as we cross the AS to DC phase boundary. 

The first fixed point FPi is reached for all points on the AS-DC surface with t greater 
than a critical value t*, and for large enough u even for t < t* . This has two relevant 
eigenvalues Ai = 4, and A2 = 3. 

The fixed point FP2 corresponds to a coexistence between the adsorbed SAW and the 
free collapsed globule phase. In Fig 18(a), the line corresponding to this fixed point is 
shown by dotted line. This point is reached if t < t*, and u is greater than, but near Uc- 

The line where the three phases meet is also an invariant manifold for the renormaliza- 
tion group flows. We find three fixed points on this line: 

1. For t<t*{= 0.34115...) the fixed point {A\ B\ S*, C*, E") = (|, i, 0.4477, 0.4528, 
0.0815) is reached. The linearized equations have three repulsive directions with 
eigenvalues Xsm = 2.2715, Ai = 3.7037 and A2 = 2.2222. The values Ai and A2 are 
the same as found for the bulk 9 point (see section 6.2). 

2. For t > t* the fixed point (i, i, 0, 0, 0.3693) is found. A gain we find three eigenvalues 
greater than one, where Xsm = 3 and the other two Ai and A2 are the same as those 
given above. 

3. For t = t* we find the symmetric "disordered and collapsed" fixed point (I;!;!?!?!)- 
This fixed point has four eigenvalues greater than one. These values are Xsmi = 
2.7620 and Xsm2 = 1.4964... and the other two are Ai and A2 as given above. 

The first two of these are tetracritical points, and the third point is an even higher 
order multicritical point ( pentacritical). 

We now discuss the behavior of phase boundaries. When t < t* = 0.34115, the critical 
line LU = uJc{u, t) which is the phase boundary between the DS and AS phases, is almost 
linear with positive slope. Beyond the multicritical point, slope of the line separating the 
AS and the DC phases rises rather sharply. In a region specified by Mc < m < (the 
value of depends on t) we have the coexistence between the adsorbed SAW and the 
collapsed globule phase. This region is shown in Fig 18 by a dotted line. For u > (t) the 
line u = LJ*{u,t) becomes almost flat. The value of Uci{t) decreases as t is increased and 
becomes equal to that of Uc t = t* = 0.34115. At t = t* the multicritical point becomes 
a symmetric "desorbed and collapsed" pentacritical point having four eigenvalues greater 
than one. 

For t > t* the critical line uj = uj*{u,t) has a different shape than for t < t* . The 
line appears to have a maximum at m < Uc- It drops rather sharply (see Fig. ITHT b) 
for t=0.5) in contrast to the case oi t < t* at the multicritical point. Furthermore, the 
hne oj = Lj*{u, t) for u > Uc separating the bulk collapsed and adsorbed phases shows the 
decreasing tendency as u is increased. The two tetracritical lines on the three-dimensional 
u — uo — t phase space meet at a pentacritical point [58]. Along one of these tetracritical 
lines, the adsorbed (swollen) polymer coexists with both the desorbed polymer and the 
desorbed globule, and the other line is the common boundary of three critical surfaces of 
a continuous transition. 
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The behavior of the special adsorption hne u = Lj*{u,t) described above can be un- 
derstood from contribution of different coexisting polymer configurations to the bulk and 
surface free energies. When both the adsorbed and desorbed phases are in swollen state, 
the adsorption line has same nature in the uj ~ u plane for all values of t, although the 
slope of the line decreases as t is increased. At t = 1, the adsorption takes place aX uj = 1 
and the adsorption line in the u — u plane has a zero slope. This is due to the fact that 
at t = 1 and u = 1 the surface is just a part of the bulk lattice. As t is increased, a; 
has to be increased to have adsorption, and since u in such a situation favors the bulk 
phase we have to increase the surface attraction to counteract this tendency. In the other 
extreme, i.e. when u >> Uc, the adsorption line has a zero slope. Here the coexisting 
polymer configurations are those given by B and E in Fig. El The free energies due to 
these two configurations balance each other at all u values and therefore the line remains 
insensitive to the value of u. It is only in the neighborhood of the special ^-point that the 
line becomes sensitive to the value of t and u. 

When t < t*, the surface layer is strongly repulsive and prohibits the occurrence of the 
E configurations in the neighborhood of the ^-point. The adsorbed state is still given by 
the configuration S, although the bulk is in the globular compact phase. Thus to balance 
the free energy u has to be increased. However, at t > t* the surface is only moderately 
repulsive and therefore at certain value of uj the polymer configuration given by E is 
formed. Thus a lower value of u is needed to balance the bulk free energy at the special 
^-point. 

A casual look at Fig. ITHT b) may give the impression of the existence of a re-entrant 
adsorbed phase as u is increased. One should, however, realize that these figures are 
merely a projection on the u-u plane of three dimensional figures in which the third 
dimension is given by t. The value of crossover exponent for the 4-simplex lattice is 
0.7481 equal to the value found for HB{2, 3). 

8.3. The 5-simplex lattice 

The grand partition function of Eq. (jH3j) for 5- simplex lattice is written in terms of 
six restricted partition functions shown in Fig 19. Out of six configurations two (A^*") 
and B^^^) represent the sums of weights of configurations of the polymer chain within one 
r-th order subgraph away from surface, and the remaining four {C^'^\ S^'^\ E^'^'^ and F^^^) 
represent the surface functions. As in the case of the 4-simplex, the recursion relations 
for A^^'^ and B^^^ do not include other variables. 

In this case, we have three phases: The adsorbed swollen (AS), the adsorbed collapsed 
(AC), and the desorbed swollen (DS). It is straight forward to write down the fixed points 
corresponding to these phases from the known fixed points for SAWs on the 4- and 5- 
simplexs. The basins of attraction of these fixed points are shown in Fig. 1201 

For the DS phase, we have A*,B* taking the value for the swollen phase of the 5- 
simplex ( section 6.2), with C*,S*,E* and F* zero. In the AS and AC phases, we have 
A*,B*, C* and E* zero, and S"*, F* has the value of corresponding the fixed points S and 
C respectively in Fig. El {^* ■, B* in the terminology of section 6.2). 

The critical surface separating the DS and AS phases is a two-dimensional surface in the 
3-dimensional parameter space (cj, m, t). All points on this surface, under renormalization, 
flow to the "symmetrical" fixed point {S* = C* = A* , E* = F* = B*), with values of 
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Figure 19. Restricted partition functions for the rth order 5-simplex. Here A^^^ and B^^^ 
represent the bulk generating functions for the polymer chain and S^'^\ C^'^\ E^'^^ and F^'^^ 
represent the surface functions. 



A*, B* same as for the DS fixed point. The linearization of recursion relations about this 
fixed point gives two eigenvalues Ai = 3.1319 and A2 = 2.5858 greater than one. The line 
Uc{u,t) is therefore a tricritical line. The crossover exponent = 0.8321. 

The AS and AC phases are separated by the critical 6'-surface uj = cudu, t) (see Fig 
20). For u; » cc;c(u, t), this surface tends to the surface u = Uc = 3.316074 in agreement 
with the critical value of u for the collapse transition in the bulk of the 4-simplex lattice. 
However, in this case, this ^-surface never meets the AS-DS phase boundary. It is shown 
in [58] that at > uj(.{u,t) the ^^-line bends and approaches very slowly the special 
adsorption line uJciu, t) as the value of u is increased. Even for u = 200 the two lines have 
not merged. The whole AS- AC boudary corresponds to the basin of attraction of a single 
fixed point, and the parameter t has no qualitative effect on the phase diagram. 

8.4. Surface adsorption of a SAW in dimensions 2 — e 

We consider the family of GM5 fractals and evaluate the value of the exponent (f) and 
examine its behavior as fo is varied from 2 to 00. The energies Eg and Et are defined as 
before. We put = for simplicity. 

A walk is called a surface walk (and assigned configuration S) when it enters through 
one corner of the surface and leaves from the other. A bulk walk represented by B has no 
step on the surface or on the bonds connecting the surface with the bulk. A walk which 
enters through one corner of the surface and ends up in the bulk is assigned configuration 
C. The generating functions for these walks can be written as 

S^'^\x,uj,t) = ^^('■)(A^,M,i?)xVi^, 

C^'\x,uj,t) = ^CM(iV,M,i?)a;Vi^. (86) 
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Figure 20. A section of the uj — u — t phase diagram for the 5-simplex at a typical value 
olt — 0.5. For other values of t, the phase boundaries shift, but the quahtative behavior 
is the same. 



Here x, as before, is the fugacity associated with each visited site of the lattice, B^'^\x) is 
the number of distinct configurations of the SAW which joins two vertices of the rth order 
of the fractal in the bulk and is the number of sites visited by SAW, M and R represent 
respectively, the number of visited sites of the lattice which lie on the surface and on a 
layer adjacent to the surface. Summations in Eqs.(85) are on the repeated indices and 
S'''^\N,M,R) and C^'^\N,M,R) represent number of configurations of respective walks 
of the r-th order gasket. The definition of parameters uj and t is same as in sec. 8.1. 
For h — 2 the following three non-trivial fixed points are found: 

1. The fixed point {B* = 0.61803,-5* = C* = 0) corresponds to the bulk state with 
V — 0.7986. For x — Xc{oj) the fixed point is reached for all u < uJc{t). Note that 
0Jc{t) is a function of t; for t — 0.5 the value of a;c(0.5)=1.1118. 

2. The fixed point {B* = C* = 0, -S* = 1) is reached for all u > ijJc{t) and x < Xc{uj). 
This represents the adsorbed state for the polymer chain with u — 1, as expected 
for a SAW adsorbed on a line. 

3. The fixed point {B* = S* = C* = 0.61803) is obtained for ui = Uc{t). The lin- 
earization of Eqs.(85) about this fixed point yields two eigenvalues greater than 
one, i.e. = 1.6709 and = 2.3819. We identify this as a tricritical point. The 

equality of B*, S* here shows that at this point the attraction at the surface, and 
repulsion at the next layer exactly compensate each other. The crossover exponent 
(j) = |°|^^ = 0.5915 At the point of adsorption transition the leading singular behav- 
ior of free-energy density is given as /(T) ~ (Tg — T)^~" Thus the "specific heat" 
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exponent a at the tricritical point is related to as ck = 2 — ^ for b — 2 the value 
of a = 0.3094. 

It is straightforward to extend this method to other members of this family. However, 
as the value of b increases the number of possible configurations of different walks increase 
rapidly. For 2 < 6 < 9 [53,54], the value of found arc 0.5915, 0.5573, 0.5305, 0.5089, 
0.4908, 0.4753, 0.4617 and 0.4497. The value of decreases as b increases and becomes 
lower than the Euclidean value of 1/2 at 6 = 6. Zivic et al [55] have used Monte Carlo 
renormalization-group (MCRG) method to obtain the value of for 2 < 6 < 100 and 
found that lower bound of Eq(84) is violated for b = 12. 

The limit b ^ oo was analyzed by Kumar et al [53] using the finite size scaling theory 
(see sec. 5) and it was found that in this limit 

<P{b) = 1 + u{b){l - D,) 

As for 6 — > oo i/(6) = |, we get 0(oo) = |. The first correction to term to finite b is 
proportional to 2 — J^, similar to that found for 7(6) in Eq(71). Since is positive 

for = |, the first correction term in Eq(86) is negative which, when multiplied by 
(1 — Db)i'{b) makes a positive contribution to (j){b). This implies that (j){b) approaches to 
I value monotonically as b is increased. 

Similar to the behavior of exponent 7 for SAWs discussed in section 5, we find 0(6) 
does not converge to the Euclidean value in the limit b ^ 00. It has, however, been noted 
in [53] that the adsorbing surface of a fractal container is similar to that of a penetrable 
surface of a regular lattice in which case 4> = 1 — 1/ should be satisfied [56, 57]. 

9. INTERACTING WALKS 

So far we were concerned with a single linear flexible polymer chain and studied its con- 
formational properties in different environments. We now show how the critical behavior 
of two interacting long flexible linear polymer chains can also be studied using a lattice 
model of two interacting walks. 

Depending on the solvent quality and the attractive interactions between intrachain 
and interchain monomers a system of two interacting long flexible polymer chains can 
acquire different configurations. The chain may be in a state of interpenetration in which 
the chains intermingle in such a way that they cannot be distinguished from each other. 
Or, the two chains may get zipped together in such a way that they lie side by side as in 
a double stranded DNA. It may also be possible, particularly at high temperatures, that 
the two chains get separated from each other without any overlap. 

By varying the temperature or tuning the interaction the system can be transformed 
from one state to another. The point at which the zipped-unzippcd transition takes place 
is a tricritical point and in its proximity a crossover is observed. In the asymptotic limit 
the mean number of monomers M in contact with each other at the tricritical point is 
assumed to behave as 



1 - 



{2u{b) - 1) log log 6 
2 log 6 



+ terms of order 



log 6 



(87) 
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where is the total number of monomers in a chain and y is the contact exponent. 

A lattice model of two interacting walks on the n-simplex lattices has been developed 
in [l^inniinni- two different situations were considered: In one, the two walks have no 
mutual exclusion, and a lattice bond can be occupied by one step of each walks. This 
model has been referred to as a model of two interacting crossed walks (TlCWs). In the 
second model, there is mutual exclusion between chains and a lattice bond can at most 
be occupied by a step of either walks. This model is called a model of two interacting 
walks (TIWs). In a general model of two interacting walks one can associate an energy 
Efy for each bond occupied by a step of both walks. The model of TICWs corresponds to 
Eb = and that of TIWs to Eb = 00. 

In the TIWs the two chains cannot be at the same site because of mutual exclusion, 
but there is lowering of energy if they occupy nearest neighbour sites. The strength of 
both the inter- and intra-chain monomer interactions depend on the solvent and chemical 
nature of the monomers. Let the two chains or walks be referred to as Pi and P2. The 
generating functions for this system can be written as 

G'(xi,X2;mi,M2,%) = E (89) 

all walks 

where Ni,xi,ui and Ri{N2,X2,U2 and R2) refer, respectively, to the number of steps, 
fugacity weight attached to each step, the Boltzmann factor associated with the attractive 
interaction between monomers and the number of pairs of nearest neighbors in the chain 
Pi{P2). R3 is the number of pairs of monomers of different walks occupying the nearest 
neighbor lattice sites and M3 denotes the Boltzmann factor associated with the attractive 
interaction between monomers of Pi and P2. 

Since the individual chain can be either swollen, collapsed or at the tricritical (9) point 
the variables {xi,Ui) and (x2,M2) can be taken to be known (see section 6.2). Therefore, 
Eq.(89) has only M3 as independent variable. With these simplifications it has become 
possible to evaluate the restricted generating functions for the n-simplex lattices. 

From the generating function of Eq.(89) one can calculate the average number of 
monomers of the two chains which are in contact (nearest neighbor) with each other 
from the relation 

{R:,} = «3^. (90) 

The RSRG transformation has been used in [60, 61] for n-simplex lattices to solve the 
model exactly and calculate the phase diagram and the value of y for the different condi- 
tions of the solvent. 

Since the topological structure of a 3-simplex lattice is such that it can not have two 
SAWs on it, the model of TIWs was solved on the 4,5 and 6 simplex lattices. These 
lattices can have two interacting SAWs, with possible self-attraction also, without the 
walks crossing each other at any lattice point. On the 3- and 4- simplex lattices TICWs 
model in a condition in which both chains are in the swollen state has been solved in [61]. 

9.1. TIWs on the 4-simplex lattice 

To describe the two walks on the 4-simplex lattice we need five restricted partition 
functions. Two of these, A^^^ and B^^^ defined in Fig 11 correspond to the walk Pi and 
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Figure 21. The state of a system of two polymer chains in a non-selective solvent on the 
4-simplex lattice in the X3{= 1/0^1X2), M3 plane. Lines SS and CC represent the tricritical 
lines of the zipped state of two chains each in the swollen state and the interpenetration 
state of the chains in a compact globule phase, respectively. Point TT, at which these 
lines meet, represents a transition point from a segregated to an interpenetrated state of 
the chains each at its 9 point. For a given value of 0:3 which corresponds to the swollen 
state of both chains, the chains are in interpenetrated state when the value of U3 is less 
than the value given by the line SS. For the value of X3 corresponding to the chains in 
their compact globule state or at their 9 points, the two chains are in the segregated state 
for all values of less than the value given by the line CC or point TT. We also show 
the X — u phase diagram of a single chain for comparison's sake. Note that the line CC 
overlaps with hne C, TT with point T. 



identical functions C^^^ and D^'^'> correspond to walk P2. The restricted partition function 
represents the configurations where walks Pi and P2 occupy neighboring sites and is 
sum of weights of configurations in which walks Pi and P2 enter and exit the r-th order 
subgraph once each. The two corner vertices of the subgraph are occupied by the walk Pi 
and the other two corner vertices by walk P2. The recursion relation for the generating 
function involves both A^'-\ S^'") and C^''\ L>(^) and is written as[60], 

£;(r+i) _ ^2(j2 ^ 2ACE{A + C) + 2E^ + QE^{B'^ + D^) + AE^{B + D). (91) 

As shown in sec. 6.2, a polymer chain on a 4-simplex lattice can be in any one of the 
states of the swollen, compact globule and the ^^-point and is described in the asymptotic 
limit by the fixed points (0.4294.., 0.04998..), (0, 22-3) and (|, \), respectively. The fixed 
points corresponding to a swollen state is reached for all values of u\ (or 1*2) < itc=3. 31607 
at xi (or X2)=Xc. The value of Xc is a function of the interaction ui (or U2). The end to 
end distance for a chain of N- monomers in this state varies as with z/=0.7294... . The 
fixed point corresponding to the compact globule state is reached for all values of ui (or 
U2) > Uc (i.e. at low temperatures) at a;i(or X2)=Xc{ui), where i=l or 2. In a compact 
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Figure 22. The phase diagram of a system of two chains in a selective solvent on the 
4-simplex lattice in the X3{= ^yxlX2), plane. Lines SC and ST respectively, represent 
the interpenetration states of two chains when one chain is in a swollen state and other 
is in compact globule phase or at its 6 point. Line CT corresponds to the configuration 
of interpenetration of the chains when one in the compact globule phase and the other is 
at its ^— point. When the value of is less than that given by corresponding lines, the 
two chains are segregated from each other. 



globule state the polymer chain has a finite density of the monomer per site when N— * oo. 
At Mi(or M2)=Mc=3. 31607 and a;i(or X2)=Xc{uc)=0.22913... the ^-point state is reached. 

In a system of two chemically different polymer chains, we may have six indepen- 
dent combinations of the individual chains which we indicate by SS,CC,TT, SC, ST 
and CT, where letters S, C and T now stand for the swollen, compact globule and 9- 
point states. The recursion relation (91) is solved using the fixed points {A*,B*) and 
(C*, D*) corresponding to the given states of the chain Pi and P2 and the starting weight 
E^^^ = (xiX2)'^ul. The solution leads to two fixed points denoted as E* and E* for each 
combination of states of the individual chains. In table I the values of these fixed points, 
relevant eigenvalues and the contact exponents y at the tricritical points are hsted [lU^. 
Figs. 19 and 20 show the state of two chains in the 0:3 (= ^/x^X2),u■i plane respectively, 
for non-selective and selective solvents. 

Lines SS and CC in Fig. 19 represent the tricritical lines of the zipped states of two 
chains each in the swollen state and the interpenetrating state of the chains each in a 
compact globule phase, respectively. Point TT, at which these lines meet, represents a 
transition point from a segregated to an interpenetrated state of the chains each at its 6 
point. For a given value of X3 which corresponds to a swollen state of both chains, the 
chains are in the interpenetrated state when the value of M3 is less than the value given 
by the line SS. For the values of 0:3, corresponding to the chains in their compact globule 
state or at their 6'-point, the two chains are in a segregated state for all values of M3 less 
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Figure 23. The phase diagram representing the configurations of a system of two chains 
on the 5-simplex lattice in the X3(= ^yxlX2), plane. Line SS, as in Fig 22, represents 
the zipped state of the chains. For all values of less than the value given by line SS, 
the two chains intermingle with each other. Line S represents the swollen state of a chain. 



than the values given by line CC or point TT. The x^u phase diagram of a single chain 
is also given in the figure for the comparison sake. 

The state of two interacting chains in a selective solvent on the 4-simplex lattice is 
shown in Fig 20. Lines SC and ST , respectively, represent the interpenetration state 
of two chains when one is in the compact globule phase or at its ^-point. Line CT 
corresponds to the configuration of the interpenetration of the chains when one is in 
the compact globule phase and the other is at its ^-point. When the value of is less 
than given by the corresponding lines, the two chains are segregated from each other. 
The individual chains configuration remains unchanged whether they are segregated or 
intermingled. 

9.2. TIW's on the 5-simplex lattice 

For the 5-simplex lattice we also need five restricted partition functions to describe the 
generating functions of two walks; two corresponding to walk Pi, two corresponding to 
walk P2 and E^'^'^ which represents the configurations in which walks Pi and P2 enter and 
exit the r-th order subgraph once each and may occupy neighboring sites. The recursion 
relation for E'^'^'> in this case is lengthy[60] and is therefore not reproduced here. 

As already shown in section 6.2 a polymer chain always remains in a swollen state for all 
values of interaction, (self-) attraction on a 5-simplex lattice. This state is characterized by 
the fixed point (0.3265..., 0.0279...). In this case we therefore have only one combination, 
i.e. SS of chains. Using these values of fixed points for single chains, the recursion 
relation for the partition function E'^'^^ of two interacting chains (see [|^) is solved and 
fixed points E*^ = 0.0279 or E* = 0.2713 are found. The fixed point E* is found for all 
values of U3 < usdxs). At M3 = u-sdxs) fixed point E* is found. 



Table 1 1 

Values of fixed points, relevant eigenvalues, and the contact exponent y at the critical point for the 4-simplex lattice. The 'a 

swollen, compact globule, and the tricritical configuration of each chain is indicated by letters S, C, and T respectively. Ai is ^ 

the largest eigenvalue of the system, and Aj refers to the two chain configuration point. ^ 

P 



State 



of individual chains 


E* 


E* 


A, 


Ai 


y 


SS 


0.04998 


0.6125 


2.6420 


2.7965 


0.9447 


CC 


0.0 


22-1/3 


2.5440 


4.0 


0.6735 


TT 


0.01484 


1/3 


2.2222 


3.7037 


0.6098 


SC 


0.0 


0.4880 


2.8520 


4.0 


0.7559 


ST 


0.02687 


0.4294 


2.8526 


3.7037 


0.6311 


TC 


0.0 


0.3680 


2.5740 


4.0 


0.6820 



fci 



Table 2 

Values of fixed points, G*, H*, I*, G}, H} and Ij. The labeling is the same as in Table 1, but for 6-simplex. 



State 








G* 


H* 




A, 


Ai 


y 


of ind. 




















chains 




















SS 


0.0175 


0.0007 


0.0007 


0.1406 


0.0147 


0.0147 


2.2054 


3.4965 


0.6318 


CC 


0.0 


0.0 


0.0 


0.0 


0.0713 


0.0713 


4.2201 


0.6318 


0.7902 


TT 


5.4x10-^ 


5.4x10-^ 


5.4x10-^ 


0.0957 


0.0535 


0.0535 


3.2225 


5.4492 


0.6902 


SC 


0.0 


0.0 


0.0 


6.75x10-5 


0.1049 


0.0711 


5.0529 


6.0 


0.9041 


ST 


0.0028 


0.0003 


0.001 


0.229 


0.1047 


0.0776 


5.2434 


5.4492 


0.9659 


CT 


0.0 


0.0 


0.0 


0.0060 


0.0714 


0.0806 


4.4019 


6.0 


0.8271 



CO 
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In Fig. 21, Ms is plotted as a function of X3 = ^/XlcX2c■ The interpenetrated state is 
found for all values of M3 which lie below lines SS at a given X3. When the values of M3 
reaches a value given by the line SS, the two chains are zipped together. Line S indicates 
the critical values of fugacity Xc and the self-attraction ui (or U2) of a polymer chain. 

9.3. TIW's on the 6-simplex lattice 

In this case one needs nine generating functions [60]; three corresponding to chain Pi, 
three corresponding to chain P2 and the remaining three denoted as G^'^\ H^^\ and J*^''^ 
represent, respectively, the sum of weights of configurations in which walks Pi and P2 
enter and exit the r-th order subgraph once each. Pi twice and P2 once, and Pi once 
and P2 twice. Since the 6-simplex lattice exhibits the collapse transition there are six 
independent combinations of single chain states, similar to the case of 4-simplex lattice. 
The fixed points corresponding to single chains found in section 6.2 have been used to find 
the fixed points of the recursion relation for the restricted partition functions for G, H, 
and / with suitable starting weights [El]- The values of fixed points , corresponding 
eigenvalues and contact exponents are given in table II. The qualitative features of the 
phase diagram found for this lattice are the same as that of the 4-simplex lattice. 

Though the model of TIW's discussed above ignores the effects of one chain on the 
critical behavior of the other chains, it provides a qualitative description of the phase 
diagram of systems of two polymer chains in a solution which may have different qualities 
for different chains and may serve as a starting point for more thorough investigation of 
segregation and entanglements in a real systems. The phase diagrams plotted here are in 
the plane xs{= y/xiX2, u^. Plot of phase diagrams in a three-dimensional Ui,U2,U3 space 
is expected to give more informations about the states the two interacting chains. 

10. CONCLUDING REMARKS 

It seems fair to say that the study of linear and branched polymers on fractals has been 
very useful in developing a better understanding of the critical behavior of polymers. 
One can verify the general qualitative features of the polymers on fractals which are very 
often similar to that in real experimental systems. One has a good deal of freedom in 
selecting the details of the fractal, and this can be used to find one that represents best 
the local geometry of the space. The exact values of the critical exponents do depend on 
the details of the fractal. But what is more important is that one can handle the different 
interactions in the problem, between different monomers, with the substrate, or with a 
different chain consistently and satisfactorily in way that allows exact calculation. 

Of course there are many unsolved problems, and possible directions for further research 
in this area. The most interesting problem would be to try to extend these exact solutions 
to some fractals with infinite ramification index. There are some studies of statistical 
physics models of interacting degrees of freedom on Sierpinski carpets, using Monte Carlo 
simulations, or approximate renormalization group using bond-moving, or other ad-hoc 
approximations. An exactly soluble case would be very instructive here. 

Determination of the exact scaling functions, not just the critical exponents, for these 
problems would be instructive: For example, the exact functional form of the scaling 
function of the probability distribution of the end-to-end distance of the SAW's on these 
fractal, or the periodic function of Fig. |H| For a numerical study of the former, see [ I64j . 
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Solving the SAW problem with quenched disorder is another interesting question. For the 
3-simplcx fractal, this corresponds to making the variables -B*^''^ random variables, and 
one has to determine the probability distribution of this variable for large r. 

We have discussed only the equilibrium properties of polymers. Of course, in many 
real systems, the time scales for cquilibriation can be very large. It is thus of interest to 
study non-equilibrium properties of statistical mechanical systems on fractals. A simple 
prototype is the study of kinetic Ising model on fractals. Closer to our interests here, one 
can study, say, the reptation motion of a polymer on the fractal substrate. This seems to 
be a rather good first model of motion of a polymer in gels. 

It is a pleasure to thank Sumedha for a careful reading of the manuscript, and her many 
constructive suggestions for an improved presentation. 
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